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Abstract 

This paper constructs the multipole expansion (in general relativity) of the 
gravitational field generated by a slowly-moving isolated source. We intro- 
duce some definitions for the source multipole moments, valid to all orders 
in a post-Newtonian expansion, and depending in a well-defined way on the 
total stress-energy pseudo-tensor of the material and gravitational fields. Pre- 
viously obtained expressions of the source multipole moments are recovered 
in the appropriate limits. The source moments parametrize the linearized ap- 
proximation of the gravitational field exterior to the source, as computed by 
means of a specific post-Minkowskian algorithm defined in a previous work. 
Since the radiative multipole moments parametrizing the radiation field far 
from the source can be obtained as non-linear functionals of the source mo- 
ments, the present paper permits relating the radiation field far from a slowly- 
moving source to the stress-energy pseudo-tensor of the source. This should 
be useful when comparing to theory the future observations of gravitational 
radiation by the LIGO and VIRGO experiments. 
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I. INTRODUCTION 



The multipole expansion is one of the most useful tools of theoretical Physics. Ex- 
tensively used in the past for dealing with wave emission and/or propagation problems 
in electromagnetism, the multipole expansion has also much contributed to gravitational 
radiation theory, starting with Einstein's pioneering work |IJ (see also 0) showing that 
gravitational radiation from a localized slowly moving source is dominantly quadrupolar, 
unlike electromagnetic radiation which is dipolar in general. 

The multipole expansion of the field generated by an isolated gravitating source can be 
said to be fully understood at the linearized approximation of general relativity. Some early 
studies extended the Einstein quadrupole formula for the energy to include higher multipole 
contributions HQ , and obtained the corresponding formulas for linear momentum and 
angular momentum J7|||. Then a number of authors PHl2l investigated the general structure 
of the infinite multipole expansion of the linearized field outside the source, and computed 
the associated fluxes of energy and momenta. It emerged from these works that the multipole 
expansion is entirely characterized by two and only two sets of multipole moments, which 
are the analogues of the electric and magnetic moments of electromagnetism and are referred 
in the gravitational case to as the mass and current moments. Notably the full expression of 
the linearized vacuum metric as parametrized by symmetric and trace-free (STF) mass and 



current multipole moments was given by Thorne [JT2 ] (we shall refer below to this expression 
as "canonical"). The use of STF harmonics is very convenient when doing computations of 



gravitational radiation [y| . 

Knowing the structure of the multipole expansion is not sufficient by itself, and one must 
relate the multipole moments to the material content of the source. In linearized gravity 



this relationship was achieved in a second series of papers [p^ -[T6[|, in which the multipole 
moments of the source were obtained as some explicit integrals extending over the stress- 
energy distribution of the matter fields, independently of any assumption concerning the 
source (it be fast moving for instance). The last of these papers, by Damour and Iyer 
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16] , put the finishing touches by expressing the multipole expansion in STF form, and by 
correcting some errors in the previous paper [ |T5|j . 

The problem is more difficult, but more interesting, in the full nonlinear theory where the 
(source) multipole moments mix with each other causing non-linear effects in the radiation 
field, for instance the so-called "tail" effect - scattering of the quadrupolar waves off the 
spacetime generated by the mass monopole of the source. Several approaches have attempted 
to define a multipole expansion consistent with the nonlinear field equations. In particular, 
when the source is stationary (existence of a time-like Killing vector), the multipole moments 
in the far field are determined by the convergent expansion of the metric at spatial infinity 
(r — > +00, t = const). Several equivalent definitions of the far-field moments are named 
after Geroch and Hansen fT7| , |T8|j , Thorne |TJ| and Simon and Beig |19 |. In the non-stationary 



case, the expansion at future null infinity (r — > +00, u = const where u is a null coordinate) 
was constructed by Bondi et al p0[ and Sachs [^TJ. By decomposing the so-called Bondi 
news function |2(| in spherical harmonics (or in STF harmonics), one obtains the multipole 
moments which are actually "measured" at infinity (radiative multipole moments). Another 
attempt at defining the multipole expansion for non-stationary sources is that of Bonnor 
and collaborators [22-^5], who combined the multipole expansion with a weak-field or post- 



Minkowskian expansion (power series in Newton's constant G). This type of approach, 
in which the multipole moments can be viewed as the source moments (distinct from the 



radiative moments), was later investigated by Thorne |12| and Blanchet and Damour |26 
(we shall refer to as paper I). 

The main result of paper I was to show that one can define starting from the Thorne 
"canonical" linearized metric an explicit algorithm for the computation of the exterior field 
up to any order in the post-Minkowskian expansion. The resulting metric represents the 
most general post-Minkowskian solution of the vacuum equations outside the source (up to 
a coordinate transformation). It is parametrized by two and only two sets of STF (source) 
multipole moments. Following paper I we shall call MPM (Multipolar-post-Minkowskian) 



this non-linear metric. It was also proved |27| that the re-expansion of the MPM metric 



at future null infinity is consistent with the corresponding expansion constructed by Bondi 



et al ||20|| . This permits defining the radiative multipole moments within the framework of 
MPM metrics. 

After dealing with the structure of the multipole expansion [^2|,|I^,^6|,^| , the next step is 
like in the linearized theory to express the multipole moments as explicit integrals extending 
over the source. At the leading order in a slow- motion expansion (speed of light c — > oo), 
the multipole moments are given by the standard moments of the Newtonian mass and 



current densities in the source (see e.g. [|12|). At the first post-Newtonian (1PN) order, 



the source moments depend on the total stress-energy distribution of the matter fields and 
the gravitational field. Because of the contribution of the gravitational field, the support 
of the total stress-energy distribution is not spatially compact, and therefore the standard 
formulas used in the linearized theory [p!^ -[T6| do not apply. If nevertheless one tries to define 
the source multipole moments using these formulas one obtains some expressions which are 
(typically) divergent due to the behaviour of the integrals at spatial infinity. This was 
the approach followed initially by Epstein and Wagoner [F2^] (see also and generalized 



formally by Thorne |T2|] to all post-Newtonian orders. The problem of divergencies becomes 
very important starting at the 1.5PN order because of the appearance of non-linear tails 
in the radiation field. Indeed one can show |3(J that ignoring the divergent terms in the 



multipole moments as defined previously is equivalent to neglecting the latter tail effects. 



Recently this problem has been solved by Will and Wiseman |3T|] who showed how to compute 
the tails within the Epstein- Wagoner- Thorne approach by means of finite expressions for the 
multipole moments. 

A satisfying derivation of the mass-type multipole moments of the source at 1PN order 
was done in |32] using a method of asymptotic matching to the interior field of the source. 



The moments were shown to be actually compact-supported, and to agree with the moments 
of Epstein and Wagoner modulo the formal discarding of infinite surface terms. The current- 
type multipole moments at 1PN order were obtained in using a similar method of 
matching. At the 1PN order there is agreement between the radiative moments (at infinity) 



and the source moments. Only at the 1.5PN order do they start differing because of the 
contribution of tails in the wave zone 



To 2PN order the multipole moments were obtained using a matching by Blanchet [30 
(we shall refer to [3(J as paper II), with a rather transparent result which seemed to be 
amenable to generalization. Namely, it was found that the expressions of the multipole 
moments depend on the stress-energy (pseudo-)tensor of the matter and gravitational fields 
in the same way as would be obtained by using incorrectly the formulas valid for compact- 
supported sources (i.e. a la Epstein- Wagoner-Thorne) , but that the multipole moments are 
endowed with a finite part operation based on analytic continuation which makes them to 
be perfectly well-defined mathematically. The latter finite part operation was found to be 
the same as used in the construction of the exterior field in paper I; it was carried in paper 
II all the way from the definition of the exterior field to the final expression of the moments. 

In this paper we shall basically generalize the result of paper II. We shall not perform 
a matching "order by order" as was done in paper II, but rather construct directly the 
multipole expansion generated by the total stress-energy pseudo-tensor of the source. This 
entails finding a formula for the multipole expansion generated by a non-compact-supported 
source. The multipole moments we shall obtain are valid for slowly-moving sources but 
to all orders in the slow-motion (post-Newtonian) parameter. Of course, the lowest-order 
post-Newtonian terms in these general expressions agree with the previous results obtained 
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32| , |33| , |30l . On the other hand, we have agreement in the limit of linearized gravity with 
the result of @. 

The multipole moments we obtain are fully consistent with the construction of MPM 
metrics in paper I as they parametrize the so-called "particular" linearized metric of paper I 
(which differs from the "canonical" metric by a gauge transformation). In a sense the present 
paper realizes a complete matching of the external field of paper I to a slowly-moving source. 
Note that the multipole moments, though allowing correctly for all the nonlinearities in the 
near-zone field, parametrize the linearized (exterior) metric which needs to be iterated in 
order to obtain the radiative moments at infinity (see paper I and [p7]). Thus the present 



results have to be combined with what we know on the relation between the source multipole 
moments and the radiative ones (see e.g. pi| , |35| and Section VI). 

Computing the moments to very high post-Newtonian order is part of a research program 
accompanying the development of the gravitational-wave detectors LIGO and VIRGO. Its 
motivation comes from the fact that in the case of the radiation emitted by compact binary 
systems an extremely precise prediction from general relativity will be necessary in order 
to extract the full potential information contained in the signal. By application of paper 
II the theoretical prediction for binary systems was obtained at the 2PN order |36| (and 
subsequently at 2.5PN order f37j). Will and Wiseman |U| obtained independently the 
result at 2PN order by application of their improved Epstein- Wagoner formalism. The work 
at the 3PN order [|38| will rely in large part on the general expressions of the multipole 
moments as derived in the present paper. 

The plan of the paper (besides this introduction and two appendices) is as follows. 
In Section II we recall some material from paper I and state our basic assumptions. In 
Section III we derive the expression of the multipole expansion of the gravitational field valid 
outside a general (slowly-moving) source. Re-writing the multipole expansion in a different 
form we obtain in Section IV the "linearized" metric (i.e. a solution of the linearized field 
equations) which is at the basis of the post-Minkowskian iteration of paper I. Decomposing 
the linearized metric into irreducible STF tensors yields in Section V the general expressions 
of the moments. The paper ends with a discussion on the link between the source moments 
and the radiative moments (Section VI). 



II. RECALLS AND BASIC ASSUMPTIONS 

The assumptions on which the present investigation is based are twofold. First we have 
the assumptions sustaining the construction of vacuum metrics by means of the Multipolar- 
Post-Minkowskian (MPM) method of paper I PB[ . The so-called MPM metrics aim at 
describing the gravitational field in the region exterior to a general isolated system. Second 
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we supplement the MPM method by other assumptions concerning the metric inside the 
isolated system. Essentially we assume that the metric is everywhere regular (smooth), 
and admits inside the system a post-Newtonian expansion which matches in the exterior 
to the vacuum MPM metric. The matching is understood in the usual sense of matching 
of asymptotic expansions. Physically the formalism is restricted to slowly-moving sources, 
whose typical internal velocities define a small post-Newtonian parameter ~ v/c. 

A. Recalls concerning the exterior field 

The Multipolar-Post-Minkowskian "exterior" metric of paper I is defined in the open 
domain Ml x M (where Ml = M? — {0}), i.e. in M A deprived from the spatial origin 
r = |x| = 0. We denote it by h^t = y/—9cxt #ext ~~ V^ 1 ' ■> where g^t is the inverse and g ext 
the determinant of the usual covariant metric, and where if" = rj^ = diag(— 1, 1, 1, 1) is 
the Minkowski metric. The MPM metric h^t is i n the form of a formal (infinite) post- 
Minkowskian expansion, 

hZ t = Gh? + G 2 h% + ... + G n /C + ... (2.1) 

(G is Newton's constant), which is such that all the coefficients of the G n, s admit a finite 
multipolar expansion in symmetric and trace-free (STF) products of unit vectors n = x/r, 
i.e. ul = STF{ni} where ni = n^n^.-.n^ (with L = i\i2...ii a multi-index with / indices; 
see |I9[ for our notation). The decomposition on the tensors h L is equivalent to the usual 
decomposition on the basis of spherical harmonics. Thus, Vn e W, 

^max 

/C(x,t) = £MMU/CM)- (2.2) 

1=0 

The expansion is considered to be finite at any given order n because in practical compu- 
tations we are doing the construction of the metric for each multipolar pieces separately. 
Note, though, that / max tends to infinity as n — >• oo. It is hoped that at the end of the 
construction it is possible to take the limit of an infinite number of multipole contributions. 
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The MPM metric satisfies the Einstein vacuum equations everywhere except at the origin 
r = [i.e. V(x, t) G Ml x M]. Using harmonic coordinates, this means 

d v hZt = , (2.3) 
□Ct = A& , (2.4) 

where the box symbol denotes the flat d'Alembertian operator □ = ri pa d p d a (with d p = 
d/dx p ), and where A^t = A^(/i ext ) is a gravitational source term, whose support extends 
out to infinity, and which encompasses all the nonlinearities, quadratic at least, of the 
field equations \K. py {h) depends on h and its space-time derivatives dh and d 2 h]. The 
relation between k. pu (h) and the Landau-Lifshitz f|| pseudo-tensor is A pu = 167iG\g\ti^/ c 4 + 
d p h pcr d a h up — h pa d 2 pu h pu . By the harmonic-coordinate condition (|2.3| ) we have, V(x, t) G 

dJCt = ■ (2.5) 

Following paper I we further assume that the metric is stationary in a neighbourhood of 
past time-like infinity, i.e. that there exists a finite instant — T in the past such that 

t < -T d/dt /COM) = • (2.6) 

This assumption permits to implement in a simple way the condition of no-incoming ra- 
diation, and, more technically, to avoid any problem of divergence at infinity of retarded 
integrals of non-linear source terms (with non-compact support). The assumption (2.6) 
could presumably be weakened in order to allow for some always radiating matter systems. 
We assume also that the metric is asymptotically Minkowskian in the past, in the sense that 

t<-T lim , xHoo /C(x) = . (2.7) 

By replacing the MPM metric (|2.1])- (|2^) into the field equations fl2.3|) - (P73D and equating 
the coefficients of the G n 's on both sides we obtain an infinite set of perturbation equations 
to be satisfied by all the h pu, s: Vn > 1, 
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QuK = , (2.8) 

Uh% = iVT • (2-9) 

For n = 1 we have Aq = 0. To any order n > 2 the source N n is known from the previous 
iterations. For instance, if A(h) = A^(/i) + 0(h 3 ) then A2 = ^2(^1)- The solution /i^t 
represents the infinite sum of the solutions of the perturbation equations (|2.8p-([2"l?D to every 
post-Minkowskian orders. Damour and Schmidt f4(| have proved that the MPM expansion 



is "reliable", in the sense that it is possible to construct smooth one-parameter families of 
solutions of the vacuum equations whose Taylor expansion when G —>■ belongs to the class 
of MPM metrics. 

The construction of the MPM metric proceeds iteratively starting from any linearized 
metric hi, solution in Ml x M of the "linearized" vacuum equations [i.e. fl2.8|) -( f2T9|) where 
n = 1 and Aq = 0] and the condition of retarded potentials. The most general linearized 
metric can be written as an explicit multipolar expansion parametrized by a set of STF 
time-varying multipole moments. This is the so-called "particular" metric of paper I we 
recall in ( ^.9| )-( P~T3| ) below. The "particular" metric differs from the "canonical" metric of 
Thorne |T2j by a gauge transformation. 

As hi is in the form of a multipole expansion it is singular at r = 0, and so is A2, 
and then successively h 2 , N n . To deal with this problem we apply when solving (2.9) 
the standard retarded integral D^ 1 on the product of N n and a factor (r/ro) s where B is 
a complex number and r a constant with the dimension of a length. The S-dependent 
retarded integral □^ 1 [(r/r ) B A r ra ] then defines a function of B which is valid (by analytic 
continuation) for all values of B in C — /Z. Near the value B = this function admits a 
Laurent expansion, and it was shown in paper I that the finite part at B = (denoted by 
FP# = o) of the Laurent expansion when B — » 0, i.e. the coefficient of the zeroth power of B, 
is a solution in Ml x M of the d'Alembertian equation with source N n . That is, we pose 

< = FP B=0 n-i \{r/r ) B N^} , (2.10) 

where d^ 1 denotes the retarded integral 



□^iV)(x,t) = ~f j^-rN&t- |x-y|/c) . (2.11) 
47T J |x — y| 



Then we have, V(x,t) e Rl x R, Uv% = jV£". 

In order to satisfy the gauge condition ( |2.8|) it is necessary to add to the particular 
solution u n a certain homogeneous solution v n whose divergence cancels out the divergence 
of u n : d v v£ v = —d v u^. Following the algorithm proposed in paper I we must decompose 
d u uf^ into four STF tensors Al, Bl, Cl and Dl and apply the equations (4.13) in paper I. 
A slightly modified algorithm, which is more convenient for our purpose, has been defined in 
the equations (2.12) of []3*5 |. In Appendix B we recall the definition of the modified algorithm 



for v n in terms of the tensors Al, Bl, Cl and Dl- Thus the solution of (|2.8| )-( ^9|) reads 
(Vn > 2) 

*C = <" + <"■ (2-12) 



When starting from the "particular" linearized metric /ipart.i [see ( ^.ll|) -( fi.l3| )1, the pre 



vious algorithm for the construction of the MPM metric generates in fact the most general 
solution of the vacuum equations under the MPM assumptions. When starting from the 
"canonical" metric /i ca n.i the MPM algorithm still generates the most general solution but 
modulo a coordinate transformation. See Theorems 4.2 and 4.5 in paper I. Crucial to the 
construction of the MPM metric is the knowledge of the general structure of the singularity 
when r — > of each of the post-Minkowskian coefficients h n , which permits to legitimate the 
application of the operator FPD^ 1 at each post-Minkowskian order. The following result 
was proved in paper I. 

Theorem. The expansion when r — > of any post-Minkowskian coefficient h n reads, V(x, t) G 
IR^x R and ViV G JN, 

h n (x,t) = J2 h L r a (\nry L F n , a , p (t) + R n , N (x,t) , (2.13) 

a<N 

where a G with — do < a < N (and ao G IN), p G IN with p < n — 1, and still / < / max [see 
(5.4) in paper I]. Both ao and / max depend on n, and tend to infinity whenn — > oo, as does 
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the maximal power of the logarithms, p max = n — 1. [The logarithms in (2.13) should really 
be ln(r/r ) but we include the constants lnr into the definition of L-F n ap .] The functions 
LF n , a ,p{t) are smooth (C°°) functions of time [starting with C°° multipole moments at the 
linearized level], and constant in the remote past: LF n , a ,p{t) = const when t < —T. They 
are given by some complicated non-linear functionals of the (source) multipole moments 
parametrizing the linear metric. The remainder i? nj jv(x, t) belongs to the so-called class of 
functions N (r N ) introduced in paper I, i.e. basically all its time derivatives are zero in the 
past (t < —T), are of class C N (M A ), and are of order 0(r N ) when r — > (we refer to paper 
I for full mathematical details). 

To prove this theorem, one assumes as an induction hypothesis that N n admits the 
same type of structure as (|2.13|) , with the only exception that p max = n — 2. Then one 



shows that applying FPD^ 1 on each of the terms composing (|2.13|) makes sense (by analytic 



continuation), and that this is a stable operation in the sense that the "solution" u n , and 
then h n = u n + v n , admits the same structure as the "source" N n , with merely an increase 
by one unit of the maximal power of the logarithms. Since fl2.13p is manifestly correct for 



hi if Pmax = 0, one concludes that (2.13) is correct (Vn) with the result that p max = n — 2 
for N n and p max = n - 1 for h n . 

B. Assumptions concerning the inner field 

We now address the problem of the gravitational field inside an actual isolated system, 
described by a stress-energy tensor T^(x, t) in some coordinate system (x, t), with spatial 
origin r = |x| = located within the system. We suppose that T^ v is stationary in the 
past, and that its support is spatially compact, T At!/ (x, t) = when r > d. The metric h^ v = 
sj—gg^ — rf v satisfies the Einstein field equations throughout M 4 , thus [using harmonic 
coordinates like in fl2.3p] 

d u W = 0, (2.14) 
UN" = l^r^ , (2.15) 



c 4 
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where we have defined the effective stress-energy (pseudo) tensor 



= \g\T^ v + -^-z=hT . (2.16) 
IQtcG 



The first term in ( |2.16 ) is the matter source term, which is of compact support. The second 



term is the non-compact-supported gravitational source term, which is the same as in 
but is computed with h instead of h ext , i.e. K^ v = k.^ u (h). From ( |2.14j ) the divergence of 



the (relaxed) Einstein equation fl2.15p yields the equation of motion of the matter fields, 

d v r^ = , (2.17) 

which is equivalent to the equation of conservation of T^ v in the covariant sense {W U T^ U = 0). 
The retarded solution of the system of equations (|2.14|) - (|2.16|) is obtained by simple inversion 
of the flat d'Alembertian operator, 

= , (2.18) 

where D^ 1 is given by fl2.11|) . Since the pseudo-tensor t^ u depends on h itself and its 
derivatives, the equation ( j2.18| ) [with the constraint equation (|2.17|) 1 should be viewed in fact 



as an integro-differential equation equivalent to ( |2.14|) -( p.l6|) and the condition of retarded 
potentials. 

We now assume that h^ v , as given by ( p. 18 ) with ( 2.17 ), satisfies certain mathematical 



properties complementing the properties of the exterior MPM metrics. Here we state our 
assumptions and comment on their adequacy in the case of the gravitational field. Based 
on these assumptions we shall construct in Section III the multipole expansion of W . 

Assumption (i) (smoothness). The field h (skipping the space-time indices) is a smooth 
function of the harmonic coordinates (x, t), namely 

h(x,t) G C°°(iR 4 ) . (2.19) 

In particular we assume no singularities in the distribution of matter fields (no point particles 
or black holes). The formalism will be a priori valid only for continuous matter distributions. 
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However, there are indications that the formalism applies also to point particles modelling 



compact objects like black holes (see p^j36|j3^] ). Simply one needs to represent the matter 
fields by delta-functions and to use a regularization in order to deal with the field near the 
point particles. 

Assumption (ii) (multipole expansion), h admits a multipolar expansion in the open domain 
exterior to the compact support of the matter source T^, in the sense that it agrees there 
with the MPM metric h ext constructed in paper I. Throughout this paper, we denote the 
multipole expansion using the script letter M.. Thus, by definition, 

M{h) = h cxt . (2.20) 

Let the exterior domain be r > TZ, where the constant TZ is such that TZ > d (with d the 
maximal radius of the compact support of T^ u ). Our assumption reads 

r>TZ M(h) = h. (2.21) 

Note that ( [2.20D simply means that we have given to /i ext the name Ai(h) (multipole ex- 
pansion of h); thus A4(h) is a solution of the vacuum equations valid everywhere except at 
the spatial origin, i.e. V(x, t) 6 ]Rl x M. By contrast (|2.21|) states that at any field point 
located in the exterior of our physical system (r > TZ > d), Ai(h) agrees numerically with h. 
Of course, h satisfies inside the system the non-vacuum field equations and therefore differs 
from the vacuum solution M.(h). For instance h is smooth all-over the system [assumption 
(i)] while Ai(h) is singular at r = 0. 

Summing up all the post-Minkowskian contributions given by ( |2.13| ) let us write the 
theorem giving the structure of the expansion when r — > of M.{h): V(x, t) G Ml x M, 
VN e N, 

M{h){x,t) = £ n L r a (\nry L F aiP (t) + R N (x,t) , (2.22) 

a<N 

where the functions L-F ap and Rn are of the type J2 n G n i 1 F n ^ v and J2 n G n Rn,N- Note that, 
in contrast with ( p,13|) , the summation integers a, p and / in (|2.22j ) take an infinite number 
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of values: a G % is not bounded from below (i.e. — oo < a < N), and p G IV, / G IV are not 
bounded from above (0 < p, I < +00). We refer to the "complete" expansion including all 
values of a G Zj in ( |2.22| ) as the "near-zone" expansion (r — > 0) of the multipole expansion, 
and denote it with an overbar, 



M(/i)(x,t) = J2h L r a (\nry L F a , p (t) . (2.23) 

This expansion should be considered in the sense of formal series, i.e. as an infinite set of 
coefficients of r a (lnr) p . 

Assumption (iii) (post-Newtonian expansion). At any fixed space-time position (x, t), h 
admits an asymptotic expansion when the speed of light c — > +00 along the basis of functions 
c _l (ln c) q where i,q G IN (h depends on c as a solution of the field equations). Thus, ViV G IN, 

/i(x, t, c) = c~' (In c) Vi,,(x, t) + T^(x, t, c) , (2.24) 

i<N 

where the coefficients <7j jg are smooth functions of (x, t) , and where the remainder T/v is of 
order 0(l/c ) when c — > 00 (non-uniformly in |x|). The corresponding expansion up to any 
order i G IN is referred to as the post-Newtonian expansion of h and denoted also with an 
overbar (this abuse of notation being justified by the matching below), 

E(x, t, c) = £ c-*(ln c) V i)? (x, t) . (2.25) 



Like ( |2.23| ) this expansion is to be viewed in the sense of formal series. Clearly the assumed 



structure Q2.24| )- (|2.25|) of the post-Newtonian expansion is consistent with ( |2.22| )-( f2.23| ). 



Indeed it is known (paper I) that the MPM exterior field, namely the multipole expansion 
A4(h), depends on the radial coordinate r only through the ratio r/c (when the multipole 
moments are considered to be independent functions of time). Thus the near-zone re- 



expansion of the multipole expansion, JA(h), can be viewed equivalently as an expansion 
when c — > 00, which should be necessarily of the type J2c~ l (hic) q [replacing r by r/c in 
(2.23)]. Our assumption ([2.25) means in fact that the multipole moments parametrizing the 



linearized metric admit themselves, when expressed in terms of the source parameters, the 
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same type of expansion ^c _:7 (lnc) s . There are many indications that the post-Newtonian 
expansion involves besides the usual powers of 1/c some (powers of) logarithms of c (see 
paper I and references therein). As usual (because all retardations r/c associated with the 
propagation of gravity at the speed of light tend to zero), the post-Newtonian expansion is 
valid only in a region surrounding the source which is of small extent as compared to one 
characteristic wavelength of the emitted radiation, i.e. h constitutes a good approximation 
to h only in the region r <C A. 

Assumption (iv) (matching). A "matching" region around the source exists, where the 
multipole expansion M.(h) and the post-Newtonian expansion h are simultaneously valid. 
In this region one expects r > 1Z and r <C A so that TZ <C A, which implies since TZ > d that 
the size of the source is d <C A or equivalently that the typical internal velocities within the 
source are v ~ dcj A C c. Existence of the matching region means therefore slowly-moving 
source. In this region we have the numerical equality [from ([2. 21 )1 



7Z<r<^\ M{h) = h. (2.26) 

We now transform (2.26) into a matching equation, i.e. an equation relating two mathe- 
matical expressions of the same nature, by replacing in the left side M.(h) by its near- zone 



expansion M.(h) as given by (|2.23| ), and in the right side h by its multipole expansion Ai(h) 
obtained from ( [2.25] ) by substituting each of the coefficients cjj i(? by their multipole expansion 
M.{a ijq ). Actually we have not defined what we mean by .M(crj >g ), as this would necessitate 
to perform a "multipolar-post-Newtonian" iteration of the vacuum equations, analogous to 
the MPM iteration of paper I. Simply we assume the existence of each A4(ai^ q ), and we 
shall determine its structure (2.30) as a consequence of the matching. Similarly to ([2.21 



the multipole expansion of h satisfies (term by term in the post-Newtonian expansion) 

r>K h = M(h). (2.27) 

Now the matching equation associated with ( p. 26 ) reads 



M(h) = M(h) . (2.28) 
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Formally this equation is true "everywhere" , in the sense that it represents an infinite set of 
functional equalities [valid V(x, t) e x M] between each of the coefficients of the series 
in both sides. Of course, the two series should be re-arranged as expansions along the same 
basis of functions, i.e. either r a (lnr) p or c~*(lnc) 9 . Note that the matching equation shows 
that the right side of (2.23) represents equivalently the "near-zone" expansion (r — > 0) of 
A4(h) and the "far-zone" expansion (r — > +00) of h. From Q2.28|) one deduces that the 
structure of the functions iF atP in ( |2.23|) in terms of the independent variable c is 

L F a>p {t, c)=J2 c~'(ln c) q L G aiPM (t) , (2.29) 

and that the structure of each multipole expansion M. (cr^g) reads 

M(a itq )(x,t) = J2n L r a (\nry L G a ^ q (t) . (2.30) 

It is important to realize that although the matching is performed in the (exterior) near- 
zone where the post-Newtonian expansion is valid (because of the small retardation r/c), 
the multipole expansion ( |2.30| ) of the post-Newtonian coefficients is valid wherever r > 1Z, 
and in particular when r — > 00. Actually (2.30) represents the far-zone expansion of cx^ g (at 
spatial infinity). Typically Ai(a iq ) blows up in the far-zone, as is clear from the positive 
powers of r in fl2.30| ). But this is not a problem, as this simply reflects the fact that the post- 
Newtonian expansion is not valid in a neighbourhood of infinity (where it would constitute 
a very poor approximation of h). 

The assumption (iv) is an application of the method of matched asymptotic expansions. 
In the present context it permits to "anchor" the multipole expansion to the field inside 
the actual source. Matched asymptotic expansions were used in general relativity originally 



for dealing with radiation reaction problems [41,42], and, within MPM expansions, in order 



to find the expression of the multipole moments with increasing post-Newtonian precision 
32] , |33|j30| , |37[] . Each time, the matching was found consistent in the sense that the satisfaction 



of the equation (2.28) determined all the desired information at the required order. 

The above assumptions (i-iv) are natural complements of the MPM framework; they 
should apply generically to the field generated by an isolated system. However let us stress 
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that we have made three physical restrictions on the system: that it be stationary in the 
remote past (before the instant — T), slowly-moving (existence of a small post-Newtonian 
parameter ~ v/c), and without singularities (no point particles or black holes). 



III. MULTIPOLE EXPANSION OF THE GRAVITATIONAL FIELD 

Our assumptions being stated, we start our investigation. Notice that the assumptions 
(i-iv) have been written for the field h, but we can readily prove with the help of the field 
equation fl2.15|) that they apply as well to the stress-energy pseudotensor r. In particular 



since Ai(A) = 1Q7tG/c 4 A4(t) [because the matter stress-energy tensor T has a compact 
support] we see that Ai(A), which is nothing but the MPM source A ext in fl2.4|) , admits 
exactly the same structure ( |2.22|) as Ai(h). Therefore we can apply on Ai(A) the finite 
part at B — of the operator □ i? 1 (r/ro) s whose definition was given in ( |2 . 1 0| ) - (|2 . 1 1| ) . We 



have seen that since Ai(A) involves an infinite sum of post-Minkowskian contributions, the 
summation in ( |2.22| ) is infinite (notably — oo < a < N). Thus we mean by FPD^X(A) the 



infinite sum of the FPD^'s acting on each separate terms composing A4(A). Recall that we 
are working in a context of approximate solutions, which are constructed only up to a given 
(though arbitrary) post-Minkowskian order, and thus involve only a finite number of separate 
contributions. Here we assume that we can consider the formal series of approximations, 



but we do not control the precise mathematical nature of this series (see, however, [|40| ). 

We can perform our derivation restricting our attention to the case where Ai(A) is not 
only constant in the past but zero in the past (t < —T) . Indeed we can check a posteriori 
that the derivation can be redone straightforwardly in the case of the constant terms in 
Ai(A) by simply using the Poisson operator A -1 instead of the retarded integral . Thus 
the constant terms can be added to the result at the end with no modification except that 
D^ 1 becomes A -1 when acting on such terms. 

From paper I the 5-dependent retarded integral □ J? 1 [(r/ro) B .M(A)] admits a Laurent 
expansion when B — > whose finite part solves the d'Alembertian equation with source 
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M(A): V(x,t) eR^xR, 

□ {FP B=0 □ /? 1 [f B .M(A)]} = M(A) . (3.1) 

From now on we pose f = r/r . Beyond the finite part all possible multiple poles ~ B~ % 
exist a priori (Vi G IV; indeed i < z max to any post-Minkowskian order n but z max — > oo 
when n — > oo). 

We want to compute the multipole expansion of the field, namely Ai(h) = 
167rG/c 4 7W(n^ 1 r). To do this we remark that because the multipole expansion satis- 
fies DM(h) = M(A) (all-over R% x Ft), we have for any B the relation U[f B M.{h)} = 
f B {M(A) + 2Br~ 1 d r M(h) + B(B + l)r~ 2 M(h)}, where d r = n^j. Applying on both sides 
of this relation the retarded integral D^ 1 and considering the finite part at B = yields 

M{h) = FP B=0 a^MiA)] 

+ FP B=0 a R 1 [f B {2Br- 1 d r M(h) + B(B + l) r - 2 M(h)}} . (3.2) 

The first term will ensure that M.(h) is a particular solution of the correct equation 
□A^(/i) = «M(A) [see (|3 ■ 1|) ] ■ This particular solution is already in the appropriate form, 



thus we concentrate our attention on the second term, which constitutes an homogeneous 
solution of the wave equation, as will follow from the fact that this second term in ( |3.9|) 
involves an explicit factor B (so it appears like a residue rather than a finite part). Because 
of this factor B we see that the contribution of any term in (|2.22| ) which is regular at r = 



will be exactly zero. This is in particular the case of the remainder in Q2.22Q which will 
give zero, and since this is true V7V G IV we conclude that A4(A) can be replaced by the 



infinite expansion A4(A) given by ( |2.23| ). Thus, writing the retarded integral in the right 
side of (ggg ) in the form ( gXLfl , 



M(h) = FP B=0 n R 1 [r B M(A)] 

- -^FP B=0 / p^[2Br- 1 d,M^+B{B + l)r- 2 Jd(h)}(y,t- |x-y|/c) , (3.3) 
47T J |x — y| 

where we denote y = y/ro. The factor B shows that the only contribution to the triple 
integral is due to a (simple) pole at B = 0, which in turn comes only from the integration 



on an infinitesimal neighbourhood of the spatial origin, |y| < e where e is an arbitrary small 
number (0 < e <C |x| say). [As we have assumed that the functions LF a , P are zero in the past 
(we add the constant parts at the end) there is no problem at the upper bound |y| — ► oo.] 
Then it is legitimate to replace the integrand in ( |3.3| ) by its Taylor expansion when |y| — > 0. 
The Taylor expansion of any K(t— |x — y|/c)/|x — y| is given by XX — ) l ULdL{K(t — r/c)/r}/l\ 
(see |$9| for the notation). In this way we obtain 

M(h) = FP B=0 D^[r B M(A)} - — £ K -^-d L \-K L {t - r/c)} , (3.4) 

where the "multipole moments" Kl{u) are given explicitely by (u = t — r/c) 

K,l{u) = -^FP B=0 / d 3 y\y\ B y L [2Br-%M(h)+B(B+ l)r~ 2 M(tij}(y,u) . (3.5) 
lb7rG J\y\<e 

Next we perform a sequence of transformations of K.l(u). By (|2.23|) , the structure of 
M.(A)(y,u) is that of a series of terms of the type X n£,|y| a (ln \y\) p LFa,p(u) ■ This shows 
that after integrating over the angles )Cl(u) is made of a series of terms of the type 
FPb=oB Jq d\y\\y\ B \y\ 2+l+a (ln\y\) p times a function of u. The latter radial integrals are 
defined by analytic continuation in B. Now the point is that the complete radial integral 
extending from to +00, i.e. Jq 00 d|y||y| B |y| 2+ ' +a (ln |y|) p , is rigorously zero by analytic 
continuation, for all values of B e C. [We repeat the reasoning already made in paper II: 
the integral can be split into the sum of two integrals, namely (d/ dB) p ^ d\y\\y\ B \y\ 2+l+a 
and (d/dB) p f^°° d\y\\y\ B \y\ 2+l+a , where A is some positive constant. When the real part 
of B is a large positive number, the first integral reads (d/dB) p {A B A 3+l+a / (B + 3 + 1 + a)}; 
and when the real part of B is a large negative number, the second integral is equal to the 
opposite — (d/dB) p {A B A 3+l+a /(B + 3 + / + a)}. These expressions represent the unique 
analytic continuations of the two integrals for all values of B except — (3 + I + a). Now 
J0" 00 ^|y| |y|' B |y| 2+i+a (l n |y|) p is defined by analytic continuation to be the sum of the analytic 
continuations of the two separate integrals, and is therefore identically zero (VI? G C).] Thus 
all the previous terms can be equivalently written as — FPb=oB / e + °° d|y| |y| B |y| 2+ ' +a (ln |y|) p . 
Furthermore, using now the presence of the explicit factor B, we see that the only contri- 
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bution to the finite part at B = comes from an arbitrary neighbourhood of the up- 
per bound |y| = +00. In this paper it is sufficient to consider as neighbourhood of in- 
finity the domain |y| > TZ, where TZ denotes the radius giving the limit of validity of 
multipole expansions. Therefore Kl(u) is given as well by a series of terms of the type 
-FP B= oBj£°° <f|y||y| B |y| 2+ '+ a (ln |y|) p , which shows that (gj5J) is in fact equivalent to 



M u ) = -77T77FPb=o / <Py\y\ B y L [2Br- l d r M{h) + B(B + l)r~ 2 M(h)}(y,u) . (3.6) 



Note the minus sign with respect to the previous expression ( |3.5|) . The next step is to employ 
our assumption (iv) of consistent matching between the post-Newtonian and multipole ex- 
pansions, according to which one can commute the order of taking the post-Newtonian and 



multipole expansions: by ( |2.28|) [or, rather, □ of (2.28)], JA(A) and M.{K) are functionally 
equal (term by term), so ( |3.6j ) can be re- written as 

-FPs=o 1 

\y\>K 

But now in the region |y| > TZ one can replace the multipole expansion by the function itself 



K l{u) = --^-FP B=0 [ d z y\y\ B y L [2Br-%M{h)+B{B + l)r- 2 M{h)]{y,u) . (3.7) 

107rG J\y\>TL 



[see ([T27D1. Thus 

— FP_B=o , 

'[y[>w 



K l{u) = --^FP B=0 / d 3 y\y\ B y L [2Br- l d r h + B(B + l) r - 2 h}(y,u) . (3.8) 
lb7rG J\y\>n 



Finally we recall that by our assumption (i) the field h and hence the function A and all 
the crag's in the post-Newtonian expansion A are regular (C°°) functions on iR 4 , and in 
particular near the origin r —>■ 0. Therefore we see that the range of integration in ( |3.8| ) 
can be harmlessly extended to the whole three-dimensional space. Indeed because of the 
factor B the contribution to the integral due to any ball of finite radius, for instance the 
ball |y| < TZ, will be zero after taking the finite part at B = 0. Thus we have 

Kl{u) = --^-FP B=0 / d z y\y\ B y L [2Br- 1 d r h + B{B + l)r- 2 h]{y )U ) . (3.9) 

Since nh = 16nG t/c 4: we can re-combine the terms in the integrand and find the equivalent 
expression 
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/ 



c 



A 



K L {u) = FP B=0 



YVl r r- 



IQ71G 



n(r B h) (y,u) . 



(3.10) 



The point about ( |3.9| ) or ( |3.10| ) is that because of the factor B the numerical values 



of the multipole moments depend on the post-Newtonian expansion of the source as inte- 
grated formally in a neighbourhood of (spatial) infinity. This is despite the fact that the 
post-Newtonian expansion is expected a priori to be valid only near the origin. As we 
have shown here, this seemingly contradictory result is possible thanks to the properties of 
analytic continuation, which permit to jump from a "near-zone" integration range in ( |3.5|) 
to the "far-zone" in (|3.6|) . However we now remark that the contributions in the complete 
multipole expansion (involving all the /Cl's) which are due to the second term in (]3.10|) 
actually sum up to give zero. Indeed we separate the d'Alembertian into a Laplacian and 
a second time-derivative, U{f B h) = A(f B h) — f B d^h/c 2 , and we integrate by parts the 
Laplacian using Ay L = 1(1 — 1)8^ _ 1 2/l-2) ) generating in this way a Kronecker symbol (the 
surface term during the integration by parts is zero by analytic continuation). Since the I 
indices are contracted with the / indices carried by the spatial gradients di = d^.-.d^ 
present in the multipole expansion [see (|3.4Q ], the latter Kronecker symbol 8^^ generates 
a Laplacian which is then equivalent (because it acts on a spherical retarded wave) to a 
second time-derivative. It is not difficult to check that the sum of all these terms with 
second-time-derivatives cancels exactly the other second-time-derivatives issued from the 
separation made above of the d'Alembertian into a Laplacian. Therefore we can re-write 
the multipole decomposition by ignoring the second term in (|3.10|) and taking into account 
only the first term which represents really the multipole moment as generated by the (post- 
Newtonian expansion of the) source. Denoting by Tii the first term in ( |3.10| ) we obtain the 
main result of this paper (restoring the space-time indices [iv): 

M{hT) = FP B=0 n R ^[r B M(An] - 7 E {-«f (t - r/c)} , (3.11) 

where the (source) multipole moments are given by the simple expression 




(3.12) 
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with T^ u the post-Newtonian expansion of t^ u . In this expression there is no explicit B- 
factor left out, so in contrast to (|3.10| ) the integration over the whole three-dimensional 
space contributes to the multipole moment, including the regions at infinity, as well as the 
intermediate regions and most importantly the near zone. As we shall see in Section IV, the 
first term in (3.11]) represents the non-linear corrections to be added in a post-Minkowskian 



expansion to the (linear-looking) multipole expansion as given by the second term. Due to 
the presence of this first term the radiative multipole moments defined in the wave zone 
will differ from the source moments (see Section VI). It can be shown that the multipole 
expansion M.(h^ v ) is actually independent of the length scale ro entering in f = r/ro and 
y = y/r [namely, the r 's cancel out between the two terms of (3.11)]. 

Let us emphasize the interesting role played by the analytic continuation all over the 
proof of (|3.12|) . Witness in particular the crucial passage from ( |3.5| ) to ( |3.6|) , which permits 
in fine to get rid of the reference to the multipole expansion in the integrand of the multipole 



moments themselves. See also the last step from (|3.10| ) to ( |3.12| ) where we discard some 
surface terms which are zero by analytic continuation, and which permits to express the 
result in terms of an integral over the sole r (without explicit reference to h). The result 



( |3.11| )-( p.l2j ) generalizes (and clarifies) the result obtained by another method in paper II 



when computing the multipole expansion of the gravitational field at the 2PN order. In the 
Appendix A we present an alternative proof of the result (|3.11|) -( [3~T2"|) based on the method 
of paper II. 



The multipole expansion ( |3.11|) - (|3.12|) has been written using non trace-free multipole 



moments. Actually it is better to re- write it using symmetric and trace-free (STF) moments, 
because the non trace-free moments are not uniquely defined (for instance Kl and Hl yield 
the same multipole expansion). Here we simply report the result of the STF multipole 
expansion (which readily follows from the equation (B.14a) in p2|): 

4Gg(-y fl fi 

1=0 

where the STF multipole moments are given by 



A(~f +00 ( _\l ( 1 ~| 

M(hn = FP B=0 n R ^ B M(An] - E -[f9 L [~^f (t - r/c) j , (3.13) 
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T^[u) = FP B=0 J d*y\y\ B y L £ dzS^z)^ (y , u + z\y\/c) (3.14) 



(jjL denotes the trace- free part of y^ [p9|). The STF moments (|3.14 ) involve an integration 



over the ^-dependent cone t = u + z\y\/c with weighting function 

(21 + IV , r 1 

Si(z) = 1 2 Zijl" ( l ~ z ) ; /_ 1 = 1 5 lim i^M*) = . (3.15) 

As a check of the result let us consider the limit of linearized gravity where we can neglect K^ v 
and .M(A^), and where t^ v reduces to the matter stress-energy tensor T^ v with compact 
support, i.e. T M1/ (x, t) = when |x| > d. In this limiting case the first term in (|3.13|) 



vanishes. Furthermore replacing by T^ u in Q3.14 ) we can remove the factor |y| s and 



the finite part at B = for compact-supported integrals. Finally, we can replace T^ u by 
T^ u within the compact support of a slowly-moving source (since d <C A). Thus we recover 
exactly the expression of the multipole moments derived for compact-support sources in the 



appendix B of [3^] and used for studying the linearized gravity in |16j (in this case the result 
is valid also for fast moving sources). 

The multipole decomposition ( |3.11| )-( p.l5| ) appears to be quite general. Physically it 



should apply to any isolated slowly-moving source without singularities. Technically it does 
not make any reference, for instance, to the post-Minkowskian (or rather MPM) expansion 
which has been invoked in order to derive it. In a sense ( |3.11| )-( [3.15| ) represents a "complete" 
matching equation (valid to any post-Newtonian and/or post-Minkowskian order), which is 
the general consequence of our matching assumption (iv). We shall leave open the possibility 
that the multipole expansion ( 3.11|) -( 3.15 ) , because of its generality, may have in fact a 



domain of validity larger than the one of MPM approximate solutions. For instance it is 
plausible that ( |3.11|) - (|3.15|) could be proved in a more general context of exact solutions 
(using perhaps an analysis similar to that of |f40|1). 



23 



IV. THE LINEARIZED MULTIPOLAR METRIC 



The most general solution of the vacuum field equations (outside the time-axis) was 
constructed within the MPM framework (see Section II). Therefore, if the present analy- 
sis makes sense, it should be possible to recast the general multipole expansion M.(h^ v ) 
as given by ( |3.13|) -( [3.14|) into a form which shows clearly that it belongs to the class of 
MPM metrics. Namely, we would like to find a certain "linearized " multipolar metric such 
that M.{h^ v ) appears to be the post-Minkowskian iteration of that metric (in the MPM 
sense). The advantage is that the multipole moments parametrizing this linearized metric, 
which constitute efficient tools in practical computations of gravitational radiation (see e.g. 
|36|-p8[), will then be obtained with full generality as computable functionals of the matter 
fields in the source. This is what we shall do in this section and the following one, where 
the linearized multipolar metric associated with M.{h^ v ) will be denoted h^rts following the 
notation of paper I. 

Let us denote the first term in ( [3.13|) by 

vT = FP B=0 n R l [f B M(k» v )} . (4.1) 

The first step is to compute the divergence d u u fMU of this term. To do this one notices first 
that M.{W) is divergenceless by (|2.14j) , and therefore that, by (|3.13|) , 

duvT = ^d u (E^^{^ r(t " r/C) l) • (42) 
If analytic continuation factors |y| B were absent in the expression of T^ 1 given by ( p,14|) , 
the right-hand-side of ( |4.2j) would be (formally) zero by virtue of d v T^ u = 0. With factors 
|y| B included, it is straightforward to find what is ( f4.2|) . In fact the computation is the same 
as the one yielding the equation (4.1) in paper II. One must evaluate the time derivative of 
J-'l' 1 using d v r^ v = 0, and perform some integrations by parts both with respect to y and 
z. We differentiate where appropriate the factor |yj . During the integrations by parts all 
the surface terms are zero by analytic continuation. Derivatives of the function 5i(z) are 
computed using 
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^[5 l+1 (z)] = -(2l + 3)z5 l (z) 



dz 



— [8 l+1 (z)\ = -{21 + 3) (2/ + !)[*,(*) - <Vi(z)] . 



(4.3a) 
(4.3b) 



As a result, we find 



'_cT 

cdu 



r L 



f -r l-i> "t 



d 



2Z + 3 \ cdu 



L ) 



(4.4) 



where ^Fi,-i> denotes the STF part of J-^-i, and where the function Q£ is given by 



g^u) = fp b=0 / d 3 y5|y| B |yr 2 ^L / ^(^(y^+zlyl/c) 



(4.5) 



With 



in hands, it is straightforward to transform 

(-)i 



into 



-t-oo C_U /-I -N 

^=^£V-M^ (t - r/c) } • 



(4.6) 



Another way of proving ( f4.5|) - ([4.6|) notices that the multipole expansion .M(A Mly ) = A^t is 
divergenceless by ( |2.5|) , and so from (|4.1| ) we have d v u^ = FPb=o □^ 1 [-Br B r _1 njA / l(A Ml )], 



where the factor B comes from the derivation of f B . Thanks to this factor the finite part is 
actually a residue, and we can perform an analysis analogous to the one going from ( |3.2|) to 
( p.9|) in the previous section. In this way we recover exactly the result (|4.5|) - (|4.6|) . 

Having obtained the divergence of u^ u in the form (|4.6|) , we proceed similarly to paper I 

a retarded solution of the 



and construct from (O) a different object v^ u which is like 
wave equation, and furthermore which is such that its divergence cancels out the divergence 
of u^ u : d u v^ u = —d u u^ u . In Appendix B we recall from paper I (and |53|]) the expression 
of v^ v in terms of STF tensors Al, Bl, Cl and Dl, and we show the equivalence with the 
different expression: 



,00 



,0; 



AG 

I 1 
AG 
~1a 



-IG° + d a (- \ c jg° a + c 2 jjg a -g, 



r 
1 • 



cjg l - -Qii 



c 



2 \r 



is: - m 



E 



^{ l -g% 
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l>2 
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(4.7a) 



(4.7b) 



l>3 



rc 
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l>2 i - 



1„„ \ _„ /l 



where the G^s are given by (4.5) (they are all evaluated at the retarded time u = t—r / c) . The 
notation is T {ij) = \ (Ty + and T| = § (rj + T/'); £^ = (0°, Q l ) means with I = 0; 
time derivatives are denoted by Q{u) = dQ(u)/du, and time antiderivatives by / Q(u) = 
J^^dvQ^v), JJ Q{u) = J^dv J Qiy). The main property of v^, i.e. d v v^ v = —d v u^ u , is 
checked directly on (|4.7| ). 

With the latter construction of v^ v we are able to define what will constitute a "lin- 
earized" approximation to the multipolar expansion J\A.{h tLl> ) given by ( |3.1^ ). Let us pose 

Aft +oo ( _\l ( 1 ^ 

Gh^ s = g ^<9 L {-JT(t - r/c)} - . (4.8) 



By the structure of /ipart.i made out of retarded solutions of the source-free wave equation, 
and by the construction of v^ u ensuring its zero divergency [see (4.2)], we have (in Ml x M) 

n/Ct.i = > ( 4 - 9 ) 
<9,/Ct.i = 0. (4.10) 

This means that (|4.8|) satisfies the linearized field equations in the exterior region, i.e. (2.8)- 
(2.9) with n — 1. By Theorem 2.1 of paper I [see also (4.7) in paper I], we know that the 
most general solution of the system of equations ( |4.9| )-( fL~T0"D can always be written as the 
sum of a "canonical" linearized metric h^ n l (introduced by Thorne ]T2J) and a linear gauge 
transformation. Therefore h^ t x [which is referred to in paper I as a "particular" metric 
but which is in fact quite general] reads 

*Cti = Ci + a*Vi + rtf - • (4-n) 

The canonical metric is parametrized by two types of (STF) multipole moments II, Jl, 
h°!.i = -4E , (4.12a) 



l>0 
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= I E {fc-i (^iL-i(«)) + y^^oL-! (^-i(«)) } . (4.12b) 
^cln.i = E pp ^4i-2(«)) + j^jd aL -2 (^e ab{l j j)bL „ 2 (u)^ j . (4.12c) 

On the other hand the vector associated with the gauge transformation depends on four 
STF moments W L , X L , Y L , Z L : 

^ = 4E^(V( W )) , (4.13a) 
4 C-V /l 

Z>0 v ' 

-pEpp (^w) + r^^-i } • (4i3b) 

All these multipole moments II, Jl, Wl, Xl, Yl, Zl will be computed in Section V. Using 
the previous definition of a linearized metric (|4.8)-( 4TT3]) , we can thus re- write our general 



multipole expansion fl3.13| ) as 

M(h^) = GU;i rtl + vT + vT (4.14) 

[where we recall that u^ u and v^ v are functionals given previously of .M(A^)]. Intuitively 
from this equation, the terms u^ u and v^ v should represent the non-linear contributions 
(of order G 2 at least) to be added to the "linearized" metric /ipart.i m order to obtain the 
"complete" (including all powers of G) multipole expansion «M(/i' u/ ). We shall prove that 
this is true, namely that by performing the non-linear iteration of /ip^rt i following exactly 
the MPM algorithm of paper I [the only difference being that we have slightly modified the 
algorithm for the computation of v^, see Appendix B], we get an infinite power series in G 
which agrees with J\A.{h lLV ) term by term in the post-Minkowskian expansion. 

Consistently with the algorithm of paper I we must consider that the first term in the 
post-Minkowskian expansion, i.e. Gh paxt _x, is purely of order G, and thus that h paxt _x itself is 
of zeroth order in G. With this convention let us show that the terms u and v in ( |4.14j ) are 
of order G 2 . Suppose they are not, so that u = Gu paurt .i + 0(G 2 ) and v = Gv paTt .i + 0(G 2 ) for 
some "linearized" coefficients Wpart.i and fpart.i- From fl4.14|) the multipole expansion of the 
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source, namely /A (A) = A(JA(h)), is given by A(G7i part .i +u + v), and therefore, inserting 
the previous assumptions for u and v, by A(G[h partl + w pa rt.i + u pa rt.i] + 0{G 2 )). Remember 
that A is quadratic in h, and pose A(h) = N 2 {h) + 0{h 3 ) [see (2.8)-(2.9)]. Thus, obviously, 
A(M{h)) = G 2 N 2 (h part ,i +M P art.i +"U P art.i) + 0(G 3 ). Now, according to (gTTJ), the finite part 
of the retarded integral of the source A(A4(h)) is u itself. Using this fact and the fact that 
the operator FPD^ 1 does not depend on G (thus it does not mix up the powers of G), we 
obtain the equation u = G 2 FPD^ 1 [A^2(/i P art.i + Wpart.i + fpart.i)] + 0(G 3 ). The right side of 
the latter equation is of order G 2 so we deduce w par t.i = 0. Then from the "linearization" of 



the formulas (|4.7|) we further deduce t> part .i = 0. Thus u and v are indeed of order 0(G 



2\ 



and thus (|4.14j) shows that M.(h^ v ) agrees with the "particular" metric at linearized order: 

M(hn = GhZ It . 1 + 0(G 2 ) . (4.15) 

Now let us denote u = G 2 u paTt .2 + 0(G 3 ) and v = G 2 v part .2 + 0(G 3 ). By the equation 
used just before (4.15), in which we can now insert M par t.i = ^part.l = 0, we have u = 
G 2 FPn« 1 [^2(W)] +0(G 3 ) thus 



= FP B=0 a^Nrih^)] . (4.16) 



Next v part . 2 is obtained from u paiTt .2 by the formulas (|4.7| ) [see also Appendix B], which are 
precisely those used in the algorithm of paper I (as redefined in ||35|| ). The definition of 
the quadratic part of the particular metric in paper I was /i part .2 = w part .2 + f pa rt.2 [indeed 
see ( p.lO|) — ( p.l2| )1, so we find from (|4.14|) that M.{h^ v ) agrees with the particular metric to 
quadratic order, 

M(hn = Gh^ tA + G 2 hZ rt . 2 + 0(G 3 ) . (4.17) 

The same reasoning is easily extended to all orders in G. 

In conclusion the "particular" metric /i part = Gh p ^ vt l + G 2 ht t ^ It 2 + 0(G 3 ) which was 
defined in paper I agrees, in the sense of powers series in G, with the general multipole 
expansion M.{h^ u ). This result is mandatory because it was shown in paper I (Theorem 
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4.2) that the general solution of the vacuum Einstein equations in Ml x M can be written as 
/ipart for some set of moments II, Jl, Wl, X l , Y l , and Z L . This is of course consistent with 
the fact that we have made no restriction when deriving the multipole expansion M.(W) 
except that it should correspond to a slowly-moving isolated system (without singularities). 
What we have gained with respect to paper I is that we understand from Section III the 
relation between M.(h^ u ) = hp^j. t [I, J, Z] and the matter distribution in the source, and 
therefore that we are able to compute the multipole moments II, Jl, Zl (given a post- 
Newtonian algorithm for the computation of t m "). 



V. THE IRREDUCIBLE MULTIPOLE MOMENTS 



By Q4.11| )-( [Ll3| ) the linearized metric h^ rt l is parametrized by six sets of irreducible 



(STF) multipole moments, with two sets II, Jl parametrizing the "canonical" linear metric 
/ican.i) and four sets Wl, Xl, Yl, Zl parametrizing a gauge transformation. We can refer to 
the moments II, Jl, Wl, X l , Y l , Z l as the source multipole moments. Of course since the 
moments Wl, Xl, Yl, Zl parametrize a gauge transformation, they do not play a physical 
role at the level of the linearized approximation. In this sense, these four moments are "less 
important" than the moments II and Jl, which constitute the "main" multipole moments 
of the source, respectively of mass-type and current-type. However, it is important to keep 
the moments Wl, • • • ,Zl as they start playing a role at the nonlinear level (at a high post- 
Newtonian approximation [[37]]). See the discussion in Section VI where we recall also that 
the six sets of moments II, Jl, Wl, ■■■,Zl are in fact equivalent physically to only two sets 
of other moments Ml and Sl- 

In the present section we use the results of Section III to compute explicitly the moments 
II, Jl, ■ ■ ■ , Zl- For this purpose it suffices to decompose into irreducible pieces the functions 
Tl* and Ql which parametrize the multipole expansions ( |3.13p and ( f4.2p . We first decompose 
the components of J 7 ^ according to 

Fl=Rl, (5.1a) 
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(o) Tl 

-l>a + Oi<i; (-) T L _i> , (5.1b) 
PI = ^U ijL + STF STF [e ain ^U ajL ^ + S Ul ^U jL ^ 

L ij 

+5ii l Saii l J~ l) U aL „ 2 + 5u ; ^ i _ 1 ( ~ 2) ^L- 2 ] + SijVz, , (5.1c) 

where the ten tensors Rl, ^T l+1 ,. . . , (~ 2 ^Ul-2, Vl are STF in all their indices. We use 
a notation standard from the works p6|Jl6|| (notably <> denotes the STF projection [j39fl). 
These ten tensors are uniquely given in terms of the J-^'s by the inverse formulas 

Rl = FT , (5.2a) 

«T L+1 = ?° L t il+1 , (5.2b) 

^T L = i ^b<L-l S H>ab , (5.2c) 

H Tl-i = f^E-i . (5-2d) 

{+2) U L+2 = TtT k+1 , (5.2e) 
21 

l~2 "L+l 



i+1) U L+ i = TT~n STF ^e H+lcd , (5.2f) 



(-D r/ _ 2(f-i)(2f-i) 

^ " (Z + 1)(2Z + 1) L-f ^-^-^ 



(5.2h) 



{ - 2) U L - 2 = , (5.2) 

^ = ^r- (5-2j) 



See for instance (5.5)-(5.8) in [16|. Next we decompose the tensors Gl according to 

Ql = Pl , (5.3a) 
Ql = {+) Q lL + e m<H ^Q L . 1>a + <W~ } ^i> , (5.3b) 



with inverse formulas 



Pl = G° l , (5.4a) 

«Q L+1 = £<;< +1 , (5.4b) 

{ ° ) QL = T ^ 1 Gb<L-i^>a b , (5.4c) 

30 



H 0m = |^i 6:i _, . (5 .4d) 

The tensors parametrizing are not independent of the tensors parametrizing . This 
is because the metric /ipart.i is divergenceless by ( f4.10|) . The (four) relations linking these 



tensors are readily obtained from (O). We have 



Pl = 


l -R L - 

c 


c 2 (2/ + l) 


(5.5a) 


M Ql = 


- (+) T L 
c 


(Z 1)^ C + l)(2i + 3) (0)f> 1 ^ 
1 J L 6c 2 /(2Z- 1)(2/+ 1) ^ c 2 (2/ + l) L ' 


(5.5b) 


(0) Ql = 


1 ( 0) f L 
c 


2 L 2c 2 (Z + 1)(2Z + 1) L ' 


(5.5c) 


"Ql = 


- Hr L 

c 




(5.5d) 



These relations permit to express the ten independent components of ft-part.i in terms of only 
six independent combinations of STF tensors. We subsitute the decompositions Q5.1J) and 
( |5.3| ) into the definition of the linearized metric fl4.7| )-( p~8| ). After some manipulations of 
STF tensors, and use of the previous relations ( |5.5|) , we arrive at a rather messy expression 
but which can be compared directly with the general decomposition of a linearized metric 
as given by the equation (2.25) in paper I. Then the six sets of multipole moments II, Jl, 
Wl, X l , Y l , Z l entering hp^ It l are obtained by applying the definitions (2.26) in paper I 
(actually our definitions differ from paper I by some constant factors). First the moments 
II and Jl are obtained as follows. In the particular cases where II has zero or one index 
(/ = 0, 1) and where Jl has one index (Z = 1), we have 

\(R + 3V) - -^f + - -JP + 4 { - } Q , (5.6a) 



c 2 c 3 c 4 c c 2 



i< = + av t ) - 4 ( " } ^ + t^ ( " 2) ^ - -JPi - H {+) Q, + ^ H Qi , (5.6b) 

J i = -- {0) T l + ^U l + 2j^Q l . (5.6c) 
c 2c z 

Then, in the generic case where II and Jl have at least two indices (/ > 2), we have 

h = + 3Vi> - -^r^n + cJ(; + 1 2 )(; + 2) '- 2 '^ • (5.7a) 

J - = + w™*' ■ < 5 - 7b > 
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Secondly, the four tensors Wl, Xl, Yl, Zl are generically obtained as 

^ = ^TT, '"'^ (» + '!)« + 2) (5 - 8a) 

^ = 2( t + lKl + 2) '^- (5 - 8b) 

y ^^TT)'^- ^ + i 2 )(/ + 2) '^- 3 ^- (5 ' 8c) 

Z L = —^Ul ■ (5.8d) 
The expressions of all these moments are obtained first by substituting into their definitions 



(|5.6|) - (|5.8|) the inverse relations (|5.2| ) and ( |5.4j ), and second by using the expressions (|3.14|) 
and ( |4.5| ) of the functions and 0%. 

We first investigate the lowest-order moments /, and Jj defined by ( |5.6[ ), which can 
respectively be called the mass monopole (or total ADM mass), mass dipole and current 
dipole of the source. As is readily checked using ( |5.5| ) (and the fact that = ( +1 ^Ui = 0), 
we have the conservation laws appropriate for gravitational monopole and dipoles, i.e. 

1 = 0, li = 0, Ji = 0. (5.9) 

For simplicity we analyse only the case of the mass monopole J; the analysis of the dipoles 
Ii and Ji is similar. The expression of / deduced from (|5.6a|) is 

Ic 2 = FP B=0 J d 3 y\y\ B £ dz{5o(r° + ¥*) - ^hy&° 1 + 

+B\y\- 2 (6 iyil T^ - c5 yjr 0i )}(y,u + z\y\/c) , (5.10) 



which can be transformed using (|4.3| ) into the simpler form 

Ic 2 = FP B=0 J d 3 y\y\ B £ dz6 {r°° - z^y^ - Bc\y\- 2 yi j¥*} (y, u + z\y\/c) . (5.11) 

The latter expression looks unfamiliar for a conserved mass, but this is simply due to the 
unusual spacelike hypersurface t — z\y\/c — u — const on which one integrates (see the 
discussion in |[L6|| ). The following technical identity is useful to transform (|5.11| ): 

^7{(^°° - ^'lyr^^OCy, ^ + = -^{lyl^Xy, ^ + . (5.12) 

32 



Multiplying this identity by \y\ B , integrating over y, and over z' from to z, and then 
multiplying by 5q and integrating over z from —1 to 1, permits to find the equivalent of 
( |5.11| ) when one uses the usual spacelike hypersurface t = u = const, 

Ic 2 = FP B=0 Jd 3 y\y\ B [r 00 - Bc\y\- 2 yjr 01 ] (y,u) . (5.13) 

The fact that / is constant is easily checked on this expression. Note that the second term 
in ( |5.13|) [whose time-derivative is associated with the flux of radiation at infinity] involves 



a factor B and therefore its value comes only from the poles of the integral at the upper 
bound |y| — > +oo. 

We now deal with the "dynamic" moments II and Jl having I > 2. In order to express 
them, it is convenient to use the following notation for combinations of components of the 
pseudo-tensor t**": 

S = 5 , (5.14a) 

cr 

S, = — , (5.14b) 
% j = r il } (5.14c) 

where r 11 = 8^7*^ . (E, E, and E^ are of zeroth order in the post-Newtonian expansion.) 
Then, by (£70, Q and ([TL4]), we find 

f f 1 \ — 4(2^ + 1) — 

h{u) = FP B = j d*y\y\ B J_ x dz |^E - c2{l + m + 3) ^iV^i 

9(91 ill ^ ) 

-8i+2VijL^ij > (y,u + z\y\/c) , (5.15) 



c 4 (/ + l)(/ + 2)(2/ + 5) 
J L (u) = e ab<il FP B =o J d 3 y|y| B J ^ dz^8iy L -i>J^b 

~ c*(i +2)(2i + 3) *w^i>«£fc}(y, u + ^iyi/ c ) • ( 5 - 16 ) 



These expressions have been derived in the non-linear theory (to all orders in the post- 
Newtonian expansion). As a check of II and Jl we can compare their expressions to the 
corresponding expressions derived by Damour and Iyer |16| in the case of the linearized 
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theory, where we can replace the pseudo-tensor T^ v by the matter stress-energy tensor T^ v 
in flat space-time (we have = T^ v inside the slowly-moving source) , and then remove the 
analytic continuation factors since T^ v is compact-supported. We find perfect agreement 
with the equations (5.33) and (5.35) in |16[]. In the non-linear theory but at the 2PN 
order, the expressions (|5.15|) -( p.l6D were already derived in paper II. To the 1PN order these 



expressions are equivalent to some different expressions derived earlier in |32],|33| (see paper 
II for the proof). 



Finally we write down the four other moments Wl, • ■ ■ , Zl- They are easily obtained as 
W L (u) = FP fc0 / i>y\y\° £ iz + 3) *"^ 

97-1-1 _l. ] 

-Si +2 yijL^ij \ (y,u + z\y\/c) , (5.17) 



2c 2 (/ + l)(/ + 2)(2/ + 5) 



r rl C 21 + 1 1 

X L (u) = FP B = j d*y\y\ B ^ dz I + + ^ y^A (y,u + z\y\/c) 



(5.18) 



Y L (u) = FP B=0 J d 3 y\y\ B £ dz { + (z + 3) ^+^' 

2(2/ + 1) — 1 

\(y,u + z\y\/c) , (5.19) 



c 2 (/ + l)(/ + 2)(2/ + 5) 



f f 1 { 2/ + 1 — 1 

Z L (u) = FP B=0 J d*y\y\ B J ^ dz I- - ^ e ab<il 5i +l y L „ l>b j: ac \ (y, u + z\y\/c) . 

(5.20) 



To Newtonian order these expressions are in agreement with the equations (4.17) in [57 



VI. DISCUSSION 



Of what use are the expressions of the STF multipole moments II, Jl and Wl, X l , 
Y L , Z L obtained in the previous equations ( |5.15|) -( [5T20"D ? From ( |4.11j )-( pl~T3| ) these moments 
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parametrize the linearized metric /ipart.i which is the "seed" of the infinite nonlinear al- 
gorithm of paper I. Thus for a specific application the expressions ( |5.15| )-( |5.20| ) have to be 
computed in a post-Newtonian expansion up to a given order and for a specific matter model 
(i.e. a specific choice of T 7 ^), and then inserted into the so-called "particular" algorithm 
of paper I for the computation of the field non-linearities (essentially this is what has been 
done for compact binary systems in |]36|-p8|). Actually the main moments to be computed 



are II and Jl because the other moments Wl, ■ ■ ■ ,Zl parametrize a gauge transformation 
and thus have no physical implications at the linearized order. In terms of a post-Newtonian 
expansion it was shown in |3(]]37j that up to the 2PN order it is sufficient to compute II 
and Jl, but that the other moments Wl, • • ■ , Zl start contributing at the 2.5PN order. 

Now it was proved in paper I that the full nonlinear metric outside an isolated system 
can always be parametrized (modulo a coordinate transformation) by only two sets of STF 
multipole moments, say Ml and Sl (different from II and Jl)- These moments parametrize 
the so-called "canonical" algorithm of paper I, defined by the same formulas as for the 
"particular" algorithm but starting with the canonical linearized metric h^ n l given by Q4.12| ) 



(but where II, Jl are replaced by Ml, Sl)- By Theorems 4.2 and 4.5 in paper I the canonical 
and particular algorithms differ from each other by a coordinate transformation. Therefore 
the multipole moments Ml and Sl are necessarily given as some functionals of the other 
moments II, Jl and Wl, • • ■ , Z L , i.e. 

M L = M L [I,J,W,X,Y,Z] , (6.1a) 
S L = Sl[I,J,W,X,Y,Z] . (6.1b) 

These functionals are quite complicated in general but can be explicitly constructed up to 
any post-Minkowskian order by implementing the coordinate transformation between the 
two harmonic coordinate systems in which the "particular" and "canonical" metrics are 
defined (see Section 4.3 in paper I). Since at the linearized level the canonical and particular 
metrics differ by a mere gauge transformation [see ( J4.ll )], we have agreement at this level 
between Ml, Sl and II, Jl- 
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M L = I L + 0(G) , (6.2a) 
S L = J L + 0(G), (6.2b) 

where 0(G) symbolizes some nonlinear (quadratic at least) products of the source moments. 
Furthermore, the result of |37|] is that in a post-Newtonian re-expansion c — > +00 we have 



M L = I L + \si L + o(\) , (6.3) 



c 5 Vc 1 - 



where 51 l is given in details by the equation (4.24) of [37 . 

On the other hand, it is known (see e.g. |jl2|) that the transverse and tracefree (TT) 
part of the spatial metric at leading order 1/R in the distance can be parametrized by still 
another double set of STF multipole moments, say Ul and Vl- These moments are called 
the radiative moments of the source, as they are the moments which would be measured at 
infinity. The radiative moments differ from the source moments because of the non-linear 
terms u^ v and v^ u in ( |4. 14 ) . Since the exterior field is entirely determined by the moments 



Ml, Sl, the radiative moments Ul, Vl are necessarily given as some (fully non-linear) 
functionals of them: 

U L = U L [M,S], (6.4a) 

V L = V L [M,S}. (6.4b) 

Ul and Vl are conveniently chosen in such a way that at the linearized order they reduce to 
the /th time derivatives of the moments Ml, Sl'- 

U L = ^f + 0(G) , (6.5a) 

^L = ^f + 0(G), (6.5b) 



where 0(G) denotes the nonlinear terms. It was shown in that the functionals 
can be constructed to all orders in the post-Minkowskian expansion by implementing the 
coordinate transformation between the harmonic coordinates and some suitable "radiative" 
coordinates in which the metric admits an expansion in powers of 1/R (without logarithms 
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of R). Furthermore, once obtained in a post-Minkowskian expansion, the functional ( p.4| ) 
can be re-expanded when c — > +00. In this limit the dominant correction is of order 1/c 3 
(or, rather, G/c 3 ) and due to the so-called tails of waves [34|| . For instance we have, in the 
quadrupole case, 



U — V\ 11 



hl| ^T + 12 



d^Mij , s „ / 1 



du A v ' Vc 5 



v) + -) , (6.6) 



where M is the ADM mass of the source (M = I), and where b is a constant time scale 
entering the relation between harmonic and radiative coordinates. The complete correction 
of order G involves other terms and can be found in |34]j35| . 

We now understand that the explicit expressions ( |5.15| )-( |O0D of the "source" multipole 



moments II, Jl and Wl,...,Zl are to be inserted into the chain of functionals ( |6.1|) and 
( |6.4| ) giving the radiative moments Ul, Vl detected far away from the source. Since by the 
present investigation the source moments have been related to the stress-energy tensor of 
the source up to any (post-Newtonian) order, and since the functionals (|6.1| ) and ( |6.4j ) can 
be "algorithmically" computed up to any nonlinear order pB|P?{| (and then be re-expanded 



when c — > 00), we conclude that the radiative moments Ul, Vl can be computed up to 
any post-Newtonian order (in principle), in terms of the source parameters. Of course the 
resulting formalism becomes extremely complicated when going to high post-Newtonian 
orders, especially when computing the source multipole moments (|5.15|) - (|5.16|) . For the 
moment it has been implemented in the case of compact binary systems up to the 3PN level 



only p-38 . 
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APPENDIX A: ALTERNATIVE PROOF OF THE RESULT 



In this Appendix we present an alternative derivation of the main result of the paper, 
which is the multipole decomposition ( 3. 11 )-( 5.12 ) or equivalently ( |3.13| )-( 3~T1 ). Though 



less elaborate than the proof presented in Section III, this alternative derivation (which 
generalizes the approach followed in paper II) permits a better understanding of why the 
multipole moments are given by ( |3. 12|) . 



We denote by A the difference between the field h, solution of the field equations Q2.14 )- 
( |2.16| ), and the finite part of the retarded integral of Ai(A) as given by the first term in 



A = h- FP B=0 n^[f B M(A)] . (Al) 



Since h is given by (|2.18| ) this difference reads 



1 (\ttC~ 1 

A = -^□flV - FP B=0 U^[f B M{A)\ . (A2) 

In the second term the finite part at B = is necessary because the multipole expansion 
Ai(A) is singular at the origin r = 0. On the other hand r in the first term of (|A2| ) is regular 
all over iR 4 , and therefore one can conveniently add into this term the finite part at B = 
without changing its numerical value (for convergent integrals the finite part simply gives 
back the value of the integral). Thus A can be equivalently re- written into the more useful 
form 

/16irG 



A = FP B=0 a- 1 



(A3) 



Under this form A appears to be the (finite part of a) retarded integral of a source with 
spatially compact support. This readily follows from our assumption that r = M.{r) when 
r > 1Z [indeed, apply □ on ( j2.21| )l, and the fact that M.[A) = 167rG/c 4 A4(r) because T has 
a compact support. Therefore the multipole expansion of A in the region r > 1Z can be 
obtained directly from the standard formula valid for sources with compact support. This 
yields immediately 
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where the multipole moments are given by 



= FP B=0 J d 3 y\y\ B y L (r - ^M(A)j . (A5) 

Now in the case of a slowly-moving source the zone of validity of the post-Newtonian expan- 
sion (or near-zone) covers the compact support of the source : d < 1Z <C A in the notation 
of Section II. Therefore both r and A4(A) in (A5) can be replaced by their post-Newtonian 



expansions, i.e. 



H L = FP B=0 I d 3 y|y|V (r- —M(A)\ , (A6) 



Finally thanks to the structure of the post-Newtonian (or near- zone) expansion A4(A) as 
given by the right side of ( |2.23| ), we see that after integration over the angles the second 
term in ( |A6| ) is a sum of terms of the type FPs=o-B J +o ° d|y| |y| B |y| 2+i+a (ln |y|) p multiplied 
by a function of time. All the latter radial integrals are zero by analytic continuation in 



B [see the discussion after (|3.5| )1, thus we conclude that the second term in flAq ) vanishes 
identically, and we recover the same result as obtained in fl3.12|) : 



n L = fp b=0 J rf 3 y |y| B y L r . (A7) 



APPENDIX B: THE HARMONICITY ALGORITHM 



The role of the tensor v^ u defined by (47/) is to cancel out the divergence of the tensor 



u iu> gj ven (41), i.e. d v {u^ v + v^ v ) = 0, while at the same time being a solution of the 
source-free wave equation, i.e. Uv^ u = 0. We first show the agreement between ( J4.7| ) and the 
expression (2.12) in which is itself a slight modification of the earlier definition (4.13) 
in paper I. The divergence d u u^ u is given by (|4.6|) in terms of the moments Q£ which are 
themselves decomposed in (|57^ ) into STF tensors Pl, ^ + 'Ql+i, ^Ql and ^'Ql-i- Posing 
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B L 
C L 



AG f-V 



Pt 



c 4 l\ 
AG (-Y -1 



c 4 l\ l + l 
AG (-)< 



(-) 



c 4 l\ 
AG (-) 1 



1 {+) Ql + 



I 



c 2 (/ + l)(2/ + l) 



c 4 Z! 



{0) Ql 



we can re-write ( |4.6|) as 

i>o v ' 

cU^ = £<9 iL (-B L 
i>o Vr 

+ ^-Cii,_i^ +e iab d aL _i 0-D 6 i_i^| . 

On the other hand the tensor t> MI/ given by (4.7) is easily transformed into 



,00 



--JA + d a (- \-cJA a + c 2 JJC a - 3B a ]) 
r \r 1 v 



-cfd + -B z 

c 



- ce iab d a ( -JD b \ - Y, 9 L-i (-A iL -i 

1 ' l>2 ^ r 



+d. 



•L-2 



1 : 

C 



3_^ 

2 " 



_ ^4j?X-2 H — CijL-2 



C 



^■ab(i^j)bL—2 ) r 



whose spatial trace is simply monopolar: 



-B 



(Bla) 
(Bib) 
(Blc) 
(Bid) 



(B2a) 



(B2b) 



(B3a) 
(B3b) 



(B3c) 



(B3d) 



We thus have agreement with the definition proposed in the equations (2.11)-(2.12) in p5 



Next we compare v^ v with the earlier definition proposed in paper I, that we denote here 
by q^ v . Using the same notation as for (|B3|) we have 



Q 



oo 



r JA + da \ r 



-cJA a + c 2 IJC a 



(B4a) 
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-.01 



~/Ci - ce iab d a [~! D b) - J2 d L-l (~ A iL-l J 7 



q ij = ~5 % 



B + d a (^B^ 



l>2 \l / \l 



+d L -2 



— 2d a L-2 [ -£ a b{iDj)bL-2 



with spatial trace 



g M = -3 

Subtracting (jB4j) from ( p3| ) we obtain 

g 00 = -39 a (±B, 



-B + <9 a 

r \r 



.00 „oo 



3 • 



v - q — — ±>i , 

rc 

^ - q v = 35^ ( -5 a ) - 6d (l ( 



which can be re-expressed in the form of the gauge transformation 



(B4b) 



(B4c) 



(B4d) 



(B5a) 
(B5b) 
(B5c) 



yen 



(B6) 



associated with the vector 



<T° = 



r 



(B7a) 
(B7b) 



Thus, had we used in ( |4.8|) the tensor q^ v instead of the tensor v i.e. had we considered 



instead of Ghp^ vt l the different linearized metric 



(B8) 



we would have obtained a dipole moment Y/ differing from the dipole Y{ as given by (|5.19 ) 
[or (|5.8c| )1 with I = 1 by the formula 
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Yl = Y t - ^ =Y i - 6 -^Q, , (B9) 

with no other modifications whatsoever. Thus the expression of the dipole moment Y- would 
not be "uniform" with the expressions of the other moments Yl for arbitrary / > 2. In this 
paper we have opted for the definition v^ u instead of q^ v for this reason, and also because 
the spatial trace v n is simpler than the corresponding q 11 [compare (|B3d|) and ( B4d| )1. 
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